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Abstract 

High-energy neutrino astronomy will be capable of observing particles at both extremely high ener- 
gies and over extremely long baselines. These features make such experiments highly sensitive to the 
effects of CPT and Lorentz violation. In this article, we review the theoretical foundation and moti- 
vation for CPT and Lorentz violating effects, and then go on to discuss the related phenomenology 
within the neutrino sector. We describe several signatures which might be used to identify the pres- 
ence of CPT or Lorentz violation in next generation neutrino telescopes and cosmic ray experiments. 
In many cases, high-energy neutrino experiments can test for CPT and Lorentz violation effects with 
much greater precision than other techniques. 
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I. INTRODUCTION 



The invariance of the product of Charge conjugation, Parity and Time reversal (CPT) is 
one of the most fundamental symmetries found in physics pj. The invariance of CPT is also 
intimately related with Lorentz invariance, another symmetry of deeply embedded within our 
current picture of the physical world. 

It is not completely clear that these symmetries will remain perfectly unbroken under all 
conditions, however. In particular, many efforts to describe the force of gravity within the 
context of a quantum theory imply the breaking of CPT and Lorentz symmetries. These effects 
are so minuscule, however, that any indication of them will likely require observations either 
at extremely high energies, or over incredibly long baselines. High energy neutrino astronomy 
provides precisely such a laboratory. 

Most varieties of particles cannot travel undisturbed over cosmological distances. Particles 
which scatter off of background radiation or other targets are of limited use to very long 
baseline experiments, particularly at very high energies. High-energy neutrinos, on the other 
hand, can travel thousands of mega-parsecs without undergoing any energy loses beyond those 
from Hubble expansion. 

Experimentally speaking, the field of high-energy neutrino astronomy is developing very 
rapidly. Current technology, such as the AMANDA II and Rice j3| telescopes at the South 
Pole, have been operating successfully for several years. The construction of next generation 
optical Cerenkov neutrino experiments such as IceCube [H and Antares is also going ahead. 
Additionally, ultra-high cosmic ray experiments such as Auger 0], EUSO [?} and OWL 8] will 
be very sensitive to the highest energy range of the cosmic neutrino spectrum. The prospects 
for other technologies which incorporate radio [§] and acoustic [HJ detectors appear to be very 
promising as well. With currently existing experiments collecting data and new experiments 
being developed and constructed, a new window into the Universe is beginning to open. 

Over very long distances, the effects of CPT and Lorentz violation can result in modifications 
to the standard neutrino oscillation phenomenology. By using high-energy neutrino telescopes 
to measure the ratios of neutrinos flavors coming from distant sources, the possible effects 
of CPT and Lorentz violation can be constrained far beyond the levels which are currently 
experimentally accessible. 

In this article, we discuss how high-energy cosmic neutrinos can be used to test for the 
violation of CPT and Lorentz symmetries. In section |H1 we review the present status of 
standard neutrino oscillation phenomenology and discuss the violation of Lorentz and CPT 
symmetries in section IIHI In section IIVI we discuss these effects in the context of neutrino 
oscillation phenomenology, and describe potential sources and experimental prospects for the 
observations of these effects in section We summarize our conclusions in section I VI II 

II. STANDARD NEUTRINO OSCILLATIONS: PRESENT STATUS 
A. Formalism 

1. Vacuum oscillations 

Normally, neutrinos are identified by their flavor (e, fi, r) rather than their mass. Neutrinos 
which have definite flavor need not be in states of definite mass, however. If this is the case, 
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we consider the flavor states, \v a ), to be a linear combination of mass states, \ui): 

W) = y^UciW%), (1) 

i 

where U a i are the components of the unitarity leptonic mixing matrix. In the Schrodinger 
representation, the time evolution of the mass eigenstates has the form 

i-r\vi(T)) = mi\ui(T)), (2) 

where r is the time in the mass frame and rrii the mass eigenvalues. Using Eq. (J2J), the 
probability of oscillation from v a to Up can be calculated to be 



P{v a ^v p ) = 5 aP - AY,^{K j U^U m U; i )^ 2 [/\m] l (L/AE)} 

j>i 

+2 £ ^(KiUpjUaJJ^ sm[Am%(L/2E)}, (3) 



j>i 



where Am|j = m 2 — mf, E is the energy of the neutrino and L is the path length. We have 
assumed the neutrino to be relativistic. If we consider three flavors of Majorana neutrinos, then 
the 3x3 mixing matrix is given by U = V ■ M, where 

C12C13 S12C13 s 13 e~ iScp 

V = [ -s l2 c 2 3 - c l2 s 2 3S 13 e iScp c 12 c 23 - s l2 s 23 s l3 e iScp s 23 c l3 

■S12-S23 — Ci 2 C 23 S\ 3 e l5cp — C12S23 — Si2C23Si3e l<5cp C23C13 

e^ 1 \ 

M= I e^ 2 . (4) 

01/ 

Here, and Sij represent the cosine and sine of the mixing angle respectively, 5 C p is a CP 
violating phase and the phases in M are the Majorana phases. For Dirac neutrinos, the situation 
is similar but the Majorana phases may be absorbed into the phases of the mass eigenstates. 
In practice, since the mixing angle #13 has been found to be small, systems involving just two 
neutrinos are often considered. Then, replacing c and h, the vacuum oscillation probability 
reduces to 

1.27Am 2i 



P{v a vp) — si n 2$ sin z 



E 



(5) 



where Am 2 is measured in eV 2 , L is measured in kilometers and E is measured in GeV. By 
convention, # 12 and Am^ are solar oscillation parameters describing v e — > whilst #23 and 
Am| 2 are atmospheric neutrino oscillation parameters, describing 



2. Oscillations in matter 

The situation becomes somewhat more complicated if neutrinos are passing through dense 
matter, for example in the Sun. In this case, the Hamiltonian in the mass basis is no longer 
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diagonal, resulting in an energy dependence in the mixing angles. For simplicity, we will consider 
a two neutrino system. Incorporating matter effects, the Hamiltonian in the mass basis is 11 



HI 



1 

2E 



m\ + A cos 2 9 A sin 9 cos 9 
A sin 9 cos 9 rrii + A sin 2 9 



(6) 



where A, embodying the matter effects, is A = 2\f2EGpN e with Gp the Fermi constant and 
iV e the electron number density. The oscillation probability therefore becomes 



where 



with 



sin 29 n 



sin 2 29 m sin 2 [fi^L] 
Am 2 sin 29 



4E Q 



M 



n 



M 



Am 2 



A 



(A V 
- — ^ - cos 29 ) + sin 2 29 
Am z 



(7) 
(8) 

(9) 



For a comprehensive overview of neutrino physics, see for example Ref. [11 



B. Where we stand - the status of neutrino oscillations 



1. Solar neutrino oscillations 



The first indirect evidence for neutrino oscillations came from those neutrinos created in the 
Sun. The standard solar model l3| describes the complex nuclear processes occurring in the 
Sun and from this, predictions of the solar neutrino flux may be extracted. Early experiments 
measuring the solar neutrino flux, such as the Chlorine 0] and Gallium [lil ]. were sensitive to 
electron neutrinos only and reported a neutrino flux significantly less than those predicted by the 
standard solar model. More recently, the SNO [l5( experiment, which is sensitive to all neutrino 
flavors, showed the total solar neutrino flux agrees well with the predictions of the standard 
solar model but with an appreciable suppression of the electron neutrino flux. This suppression 
is definitive evidence of neutrino oscillations from u e — > v^ T . The KamLAND experiment 
which detects electron anti-neutrinos from nuclear reactors, reported a significant suppression 
in event rates thus corroborating the SNO results. Using data from these experiments, the 



solar neutrino oscillation parameters, Am^ 
with high precision. These values are shown in table 



Amg and sin 2 9 
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sin #0, have been measured 



2. Atmospheric neutrino oscillations 

Whilst the first indirect evidence for neutrino oscillations came from solar neutrinos, the 
first direct evidence came from atmospheric neutrinos. Atmospheric neutrinos are created by 
interactions of cosmic rays with atmospheric atomic nuclei. Pions, created in this interaction 
decay by 7r~ — > + ]7 M ^ + + u e + e~ (with an analogous decay chain for 7r + ), thus 
indicating that the muon neutrino flux should be roughly twice that of electron neutrinos. 
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Parameter 


Value 


2/3 neutrino analysis data 




8.3 x 10" 5 eV 2 


2 


solar + KL(766.3 TY) 


A 2 


2.3 x 10 ev 


2 


SK + K2K 


sin 2 00 


0.27 


2 


solar + KL(766.3 TY) 


sin 2 otm 


0.50 


2 


SK + K2K 


sin 2 6»i 3 


< 0.05 


3 


solar+KL(766.3 TY) 
+ CHOOZ + atm 



TABLE I: Current status of neutrino oscillation parameters found from the Super-Kamiokande (SK), 
K2K, KamLAND, (KL) and CHOOZ solar and atmospheric neutrino experiments. These values have 
been taken from Ref. [2jJ. 



3 and Soudan jl^] experiments reported a significant 



12(1 l2lL |22f which suggested oscillations from to v T . 
however. In 1998, this changed when the Super- 



However, the Kamiokande [171 ] . 1MB 
deficit in the expected z/ M : v e ratio 
These results were somewhat inconclusive 
Kamiokande collaboration showed the muon neutrino flux had a zenith angle dependence, which 
implied a dependence upon path length j23|. This was the first direct evidence that neutrinos 
oscillate from one flavor to another. A summary of the present values of the atmospheric 
neutrino oscillation parameters is given in table ED 



3. Three neutrino oscillations 



The standard analysis of solar and atmospheric neutrino flux is done within the two neutrino 
approximation. However, in order to place values on Am 2 x and sin 2 # 31 , the data must be 
analyzed taking all three neutrino flavors into account. Table H] again shows the present status. 



4. The LSND result 



It would seem from the discussion above that, short of higher precision measurements of 
the oscillation parameters being performed, the phenomenology of neutrino oscillations is well 
understood. However, this is not entirely the case. The results of the LSND experiment |24j . 
which produces a beam of v e , and T7 M and then searches for the appearance of V e that have 
oscillated from 77^, cannot be reconciled with the solar and atmospheric neutrino data. The 
LSND results imply a mass difference which lies in the range of 0.2 < Am 2 LSND < 10 eV 2 
[25l |. If this result were corroborated by the miniBoone experiment [26], then it would provide 
indications of new physics. In particular, in order to combine in a compatible way the LSND 
result with the atmospheric and solar oscillation data, one would have to invoke oscillations 
into sterile neutrinos or break CPT invariance. This latter possibility is discussed below. 



III. VIOLATING LORENTZ AND CPT INVARIANCE 



The issues of CPT invariance violation (CPTV) and Lorentz invariance violation (LV) are 
intimately related. The CPT theorem is a fundamental ingredient of quantum field theory 
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ensuring that quantities appearing in the theories, such as the Hamiltonian and Lagrangian 
density, are invariant under the combined operations of charge conjugation (C), parity reflection 
(P) and time reversal (T). The CPT theorem holds in flat space-times provided the theory obeys 

• locality, 

• unitarity, 

• Lorentz invariance. 

Deviation from any one of these requirements leads to CPTV. It has also been shown recently 
that CPTV leads to LV Q. 

A. Violations of Lorentz invariance 

The breaking of Lorentz symmetry may arise as a consequence of quantum gravity from non- 
trivial effects at the Planck scale, such as the existence of a fundamental length scale. Naively, 
one would expect this length to be the same in all reference frames due to its fundamental 
nature. This invariance between frames is in direct disagreement with special relativity, the 
Lorentz transformation predicting length contraction. Thus, there are three possibilities for the 
fate of LV. Firstly, Lorentz invariance may hold at the Planck scale, so that the flat space-time 
picture is valid. In this case, it is our naive thought experiment which is wrong. Secondly, 
Lorentz symmetry may be broken at the Planck scale suggesting a class of preferred inertial 
frames. Very often, this preferred frame is assumed to be the related to the cosmic microwave 
background. In this case, the dispersion relation for energy and momentum is modified and 
depends upon the quantum gravity environment. Thirdly, Lorentz symmetry may be deformed 
[29]. Again, this leads to a modified dispersion relation but with the Lorentz transformations 
now containing a second observer independent scale. For example, the Planck length could be 
an independent scale in addition to the speed of light. 

1. Lorentz invariance violation in string theory 

String theory approaches the quantum gravity problem from a particle physics perspec- 
tive. Although string theory does not quantize space-time since it describes the background 
space-time entirely classically, the issues related to this are still far from being resolved. At 
this moment in time, since the background is classical, there are no indications that Lorentz 
invariance is broken within string theory. Of course, if it is found that the background space- 
time needs to be quantized, then this could lead to Lorentz invariance violating string theories. 
Having said that, two particular theories which are considered to be low energy limits of string 
theory, namely flat non-commutative space-times [3(J and the Standard Model Extensions 
(SMEs) j3l|, indicate the presence of LV. The non-commutative space-time approach assumes 
that space-time coordinates do not commute, leading to the breaking of Lorentz symmetry 
and various forms of the dispersion relation [i^, H^. The SMEs extend the Standard Model 
Lagrangian to include all LV and CPTV operators which are of dimension 4 or less, in order to 
be renormalizable. Again, this phenomenological model modifies the dispersion relation. 
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2. Lorentz invariance violation in loop quantum gravity 



In contrast to string theory, loop quantum gravity (also known as canonical quantum gravity) 
approaches the quantum gravity problem from a general relativity perspective. The theory is 
fully background independent, as is general relativity, which leads to the prediction of discrete 
space-times j34|. Initially, it was thought that loop quantum gravity would preserve Lorentz 
symmetry, however it is now thought that this theory could break j3^, 3^, 3?J or deform j^l, l39l | 
Lorentz symmetry, thus leading to modified dispersion relations. The issue of whether Lorentz 
symmetry is preserved, broken or deformed is still unresolved in loop quantum gravity. 

B. Violations of CPT 

The breaking of CPT invariance may occur independently of LV effects from the loss of 
unitarity leading to quantum decoherence. Again, this loss of unitarity could occur because of 
the discrete and topologically non-trivial nature of space-time resulting in the vacuum creation 
of quantum black holes with event horizons having radii of order the Planck length. This 
continuous creation and evaporation of these quantum singularities results in space-time having 
a foamy nature. When particles pass by these quantum black-holes, some of the particles' 
quantum numbers could be captured by the space-time fluctuations. With the evaporation of 
the black holes, the information captured would be lost to the vacuum, inaccessible to low energy 
experiments. This loss of information would imply that initially pure quantum states may evolve 
into mixed quantum states; a process forbidden within standard quantum mechanics. Since we 
are now dealing with mixed quantum states, we use the density matrix formalism. Including 
decoherence linearly, the time evolution of the density matrix, p, is modified: 

f t =- i [H,p]+V[p], (10) 

where H is the Hamiltonian of the system interacting with the environment through the op- 
erators Dj,D^. We would expect the the CPTV term, T>[p\, to take the Lindblad form [4Qj] 

V[p] = J2 fa n :, n A ~ 2D i D j) ( n ) 

j 

where {...} represents an anti-commutator. From a physical point of view, we require energy 
conservation and monotonic increase in the von-Neumann entropy. In this case, we find we have 
to specify the operators, D, to be self-adjoint and that they commute with the Hamiltonian. 
We therefore find 

V[p\ = Y J \D 3 AD v p}}. (12) 

3 

From a quantum gravity perspective, we would expect the operators, D, to be proportional to 
the inverse of the Planck mass and thus T>[p] oc M~ 2 . 

1. Quantum decoherence in string theory 

The evolution of mixed states into pure states, as described above, results in problems 
with defining an 5* matrix. Since string theory relies on the defining of S matrices, quantum 
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decoherence is generally not expected within string theory. However, one class of string theories, 
namely non-critical string theories, may allow decoherence. This theory may be viewed as a type 
of non-equilibrium string theory with the so called critical strings corresponding to equilibrium 
points within the theory (for more details see, for example, Ref. In this case, we find an 

analogous expression for the time evolution of the string matter density matrix as with Eq. (|12j) 



2. Quantum decoherence in loop quantum gravity 

Whilst loop quantum gravity implies that space-time is discrete, there is no a priori reason 
to expect quantum decoherence. However, there have been proposals jZiJ suggesting that 
the discreteness of space-time may induce decoherence having the Lindblad form outlined in 
Eq. (|10[). However, it seems that there is still much work needed in order to clarify this. 



3. Cosmological decoherence 

In addition to quantum decoherence induced by space-time foam effects, it may be that 
decoherence arises from cosmological considerations. It is now established that the universe has 
entered a period of acceleration ji^, driven by some exotic dark energy. If this expansion 
continues, the universe will evolve into a de-Sitter universe, expanding at an exponential rate. 
This would imply the existence of a cosmological horizon. This situation can be considered 
in the same way as that of the space-time foam situation except we, as the observers, now 
inhabit the space within the horizon instead of outside. The existence of this horizon would 
again lead to the inability to define S matrices leading to decoherence. It has been argued 
in non-critical string theory j4?| that this cosmological decoherence may be intimately linked 
with quantum gravity. Considering a two level neutrino system, the cosmological decoherence 
parameter, 'jcosmo, is related to the cosmological constant, A, the weak string coupling, g s , the 
difference of the squares of the mass eigenstates, Am 2 , the energy of the neutrino, E, and the 
string mass scale, M s : 

Icosmo ~ E 2 M ' 

IV. NEUTRINO OSCILLATION PHENOMENOLOGY WITH CPTV AND LV 
A. Neutrino oscillations and LV effects 

1. Modified dispersion relations 

As discussed above, if we allow LV, then this leads to modified dispersion relations (MDR). 
From the discussions above, we find it useful to parameterize the MDR's, to leading order in 
the Planck energy, Ep, as 

E 2 = p 2 + m 2 + r] P 2 (Jpj (14) 
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where E is the energy of the neutrino, p is its momentum, m is the mass eigenstate and rj and 
a are LV parameters. Assuming that that the parameter rj is not universal and depends upon 
the mass eigenstate, we may write the two neutrino Hamiltonian in the mass basis as 

ff -("V (i5) 

where we have neglected the kinetic term and identified p with E since the mass eigenstates 
and the LV are terms are much smaller than the momentum of the neutrinos. For simplicity, 
we have absorbed the Planck energy into rj. This leads to the neutrino oscillation probability, 



P[v a —* v p\ = sin 2 29 sin 2 



Am 2 L Ar/E n L 
AE + 4 



(16) 



where n = a + 1 and Arj is the difference between the two values of rj. If there are no LV 
effects, we recover the standard neutrino oscillation probability (jSJ) (replacing c and K). We 
also assumed that the LV parameter, 77, had a dependence on the mass eigenstate. If this is not 
the case, then the neutrino probability remains invariant even if Lorentz invariance is violated 
or deformed. Assuming that these effects take place in atmospheric neutrino oscillations, the 
LV effects become significant when 

L27 Am 2 L ^ i x w 27 Ar]E 2 L (17) 
E 

where we have set a = 1 for simplicity. We therefore find 

A " ~ Mr ~ irS ° eV " (18) 

using the value for Am 2 from table H] and E = 1 GeV, the peak in the atmospheric neutrino 
flux. 



2. Neutrino oscillations and the SME 

As we outlined above, the Standard Model Extension (SME) can be considered as a low 
energy phenomenological model of string theory. The effective SME Hamiltonian describing 
flavor neutrino propagation, to first order, is |49j 

K f p f = \Jt\8ap + ^[m 2 + - (c L r P ,p v )] a p (19) 

z\ p 1 

where in is related to the standard neutrino mass, a, j3 are flavor indices and Ol, cl violate 
Lorentz invariance and CPT invariance. One of the main differences between this model and 
the LV effects described in the last section is that the Hamiltonian need not be diagonal. In 
Ref. [ill, various assumptions are made in order to simplify the model. Here, instead, we will 
adopt a general off-diagonal formalism. 

In the two neutrino case, we will assume a Hamiltonian in the mass basis of the form 



H.n=\ S °i ^J. ( 2,» 
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where a\ and 02 are real, off- diagonal, LV parameters (the a's here are independent of the a 
in Eq. (|19))). In order to calculate the probability, we will use the density matrix formalism. 
Writing the Hamiltonian in terms of the Pauli matices gives 





Am 2 
4B 


Am 2 


— CL 2 

a 1 





«2 


— CLi 






^ = - 2 W Ql (21) 
\ a 2 -a a / 

where we have omitted the zeroth components for simplicity as they are all identically zero. 
This matrix has eigenvalues, Aj, given by {±zO, 0} where Vl = \J u 2 + a\ + a 2 and = Am 2 /4£ 
and the matrix in Eq. (}2T]) is diagonalized by the unitary matrix 



1 



t/ = — _ J ioa 2 + ia\VL ua 2 — ia\VL —a 2 \/2a\ + a 2 , ] . (22) 

V'2a 1 + 2a 2 fi^ fl2 + fl 2 fl2 + fl 2 + a 2 



The components of the density matrix are given by 

Pi (L) = ^^e A ^f/r fc V fe ( ), (23) 

where Uy are components of the matrix in equation ()22|) and p(0) is the density matrix initially. 
Assuming we have a muon neutrino which oscillates into a tau neutrino, the probability of 
oscillation is given by 

P = Tr[p,(L)p T (0)] (24) 

with 

,„ N / cos 2 9 sin 9 cos 9 \ / sin 2 — sin#cos#\ 

/M°) = • fl fl • 2 a ' Pr(0)= 2 . (25) 

\ sin 9 cos 9 sin / \ — sin 9 cos cos 9 J 



Thus the probability of oscillation is 



os 2 2*(l-^-^cos(2^)) 
+ sin 2 2£ (l-^-i^±^cos(2fiL)) 
-^in4^-^isin 2 (fiL))' 



(26) 



with a = a\ + zg^. In an analogous way to the diagonal case, the LV parameter, a, may have 
an explicit dependence on the neutrino energy. To examine this energy dependence, we let 
a — > aE n where n is an extra parameter of the theory. It is also particularly interesting to note 
that if we let the quantities in Eq. (J2*U|) go to 



m\ 


— > m\ 






m\ 


— > m 2 




A 


a 1 


~^ 2E 


a 2 


= 0, 



,2 



sin 9 cos 9, 

(27) 
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then we recover the Lorentz invariant matter effects situation as described in section |H] 

If we wish to include these off-diagonal LV effects in the case for three neutrinos, the situation 
becomes very difficult as, we have three mixing angles, three mass differences and three LV 
parameters. In order to examine how the LV effects manifest themselves, we consider only 
a first order approximation in the LV parameters. In the standard oscillation case, the time 
evolution of the density matrix is given by 

where B is the matrix representing the Hamiltonian in the Pauli basis. Perturbing the density 
matrix and the matrix B: 

p -> p + Spx 

B -> B + 5C, (29) 

where the 5 quantities contain the LV effects. Substituting into (|28|) and equating coeffi- 
cients gives 

8px = B5 Pl + 5C Po . (30) 

Defining the vectors x, y as 

Po = Ux, 5pt = Uy, (31) 

where the components of the p vectors are the components of the density matrix and U is the 
unitary matrix diagonalizing B, we may rewrite Eq. (fHUJ) as 

y - U~ x BUy = IJ- 1 5CUk. (32) 

Since we know U, B, 5C and can evaluate x, then solving this equation will give us the 
perturbation to the density matrix from which we may calculate oscillation probabilities. In 
reality, this calculation still results in complicated expressions for the probabilities. However, 
the expressions are greatly simplified if we assume very long path lengths. This is entirely 
reasonable since we need only consider the three neutrino system when considering neutrinos 
from astrophysical sources. Using the values of the mixing parameters from table HI assuming a 
normal mass hierarchy and that the LV parameter, a, is real, we find the oscillation probabilities 
to be 



p 


[V e - 




= 0.564- 4.39 x 10 9n+12 aE" +1 , 


p 


~y e - 


+ "*J 


= 0.264 + 1.54 x 10 9n+12 aE n+1 , 


p 




-> V T \ 


= 0.180 + 2.93 x 10 9n+12 a£ n+1 , 


P[ 


^ - 




= 0.365 - 1.30 x 10 9n+11 aE" +1 , 


p 




-> U T ] 


= 0.367- 1.16 x 10 9n+12 a£" +1 , 


p 


v T - 


-> V T ] 


= 0.449 - 1.56 x 10 9n+12 aE n+1 , 



(33) 

where we have written out in full the explicit dependence of the LV parameter on the neutrino 
energy and replaced c and H. These results are only valid near the threshold of LV effects setting 
in. Using these probabilities, it is possible to find expressions describing the flux of neutrinos 
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originating in astrophysical sources. If we assume that only electron and muon neutrinos are 
created, we parameterize the initial flux as 



= (1 (34) 

where e G [0, 1] and $>tot is the total flux. In terms of the neutrino probabilities of Eq. (J33)) . the 
neutrino flavor composition at the detector is given by 

R Ve = (P[u e -> v e ]$ Ve + -> v e \<S> Vfi 

tot; 
tot; 

i2„ T = (P[i/ e -> i/ T ]$„ e + ^1$^ 

+ P[v T -> z/ T ]$, T )/$ tot (35) 



and so we find 



i2„ e = 0.264 + 0.30(k - aE n+1 [0.593e - 0.154] x 10 9n+13 , 
R Vfl = 0.365 -0.101e + aE n+1 [0. 167 e - 0.013] x 10 9n+13 , 

= 0.367 - 0.187£ + aE n+1 [0A09e - 0.116] x 10 9n+13 . (36) 

Again, this result is only valid near the threshold of LV effects. We have also calculated the 
ratios in the large a limit and find: R Ue : R Vjl : R Vt ~ 0.42 : 0.57 : 0.013 for e — 1/3 and 
R Ue : R Vfi : R Vt ~ 0.70 : 0.27 : 0.027 for e = 1. Our numerical calculations indicate that the 
transition between standard oscillation phenomenology and the phenomenology of the large a 
limit takes place suddenly. If we assume what is perhaps the most natural choice of a = M^ 1 
(n = 1) or a — Mp L 2 (n — 2), the thresholds for these effects take place at ~ 1 TeV and ~ 10 
PeV, respectively. Approximate numerical results are illustrated in figure [T] 



B. Quantum decoherence and neutrino oscillations 

We now turn to the violation of CPT without explicit Lorentz violation and consider how 
quantum decoherence affects neutrino oscillations. As we discussed in the previous section, 
quantum decoherence causes the time evolution of the density matrix to be altered: 

d f t =-i[H,p\+5$p (37) 

where arises due to the loss of coherence. Considering two neutrinos and expressing this 
equation in the Pauli matrices basis, we find the time evolution of the density matrix to take 
the form 

^ = (V + K v ) Pv (38) 
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FIG. 1: The estimated effects of Lorentz violation on the ratio of neutrino flavors observed after 
propagation. The top frames correspond to neutrinos produced through pion decay, while the lower 
frames correspond to (anti-)neutrinos produced through the decay of neutrons. The left and right 
frames display the behavior of models with Lorentz violating parameters proportional to E 2 and E 3 , 
respectively. 



where h represents standard oscillations and h' the decoherence effects. Following Ref. [50], we 
parameterize h' as 

/o \ 

h' = -2 ° a b d . (39) 
b a (3 y ' 

\0 d p 5 j 

The first row and column contain only zeros in order to conserve probability and obey the 
second law of thermodynamics. We may further constrain the theory if we assume that energy 
is conserved within the neutrino system. In order for this to be the case, the parameters d, (3, 5 
must all be identically zero. However, it is not clear what assumptions we should make from a 
quantum gravity perspective and so we will also examine the possibilities that these parameters 
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have non-zero values. Using Eqs. (|S5|l and we obtain the equations, 



Po 
Pi 

f>2 
P3 



(. Am 2 \ 
-2ap l -216- I p 2 - 2dp 3 , 

/, Am 2 \ 
-2 ( 6 +^J Pi-2ap 2 -2/5p 3 , 

-2dpi - 2/3p 2 - 25p 3 . 



(40) 



In order to find the oscillation probability, we solve these equations with suitable initial condi- 
tions. There does not exist, however, a simple closed form for the solution to these equations 
and so, for illustrative purposes, we consider two limiting cases. The first is that the parameters 
d and f3 are zero and the second is that all the decoherence parameters are zero except d. In 
the first case, the oscillation probability takes the form 



PK^M = -{cos 2 20(1 -e 



-2SL\ 



+ sin 2 20 



1 - e - {a+a)L cos 2L 



Am 2N 



AE 



4 V ; 



i \ -] 

2 



(41) 



If we assume complete positivity |51| . necessary for the Lindblad form, (as described in section 
IIII|) with energy conservation, we find the simplest possible extension to standard neutrino 
oscillations which includes dehoherence. In this case a = a with all other parameters equal 
to zero. Considering the second approximation, with non-zero d only, gives the probability of 
oscillation to be 



P\v n 



. cos 2 29 
- \ 

+ sin 2 29 



uj 2 d 
+ 



n2 ' fi2 cos(20 d L) 



+ sin 40 



d 

d 2 — uj 2 
1 + Q2 cos(20 d L) 

— sin(20,L) 



(42) 



where Vtd = \Jw 2 + d 2 and uj = Am 2 /4£'. Note the similarity to the probability given in 
Eq. (HSJ). 

So far, we have said nothing about the form of the decoherence parameters. It is clear 
that the decoherence parameters must have dimensions of energy but it is also possible that 
they have an explicit dependence on the neutrino energy. In the literature [52J, |53|, |54j , three 
models have received significant attention, specifically, with energy dependences E°, E -1 and 
E 2 . From a quantum gravity perspective, an energy dependence of E 2 is particularly interesting 

mm. 



In a similar way to the LV case, it is also worthwhile studying the three neutrino case. 
Again, the situation is made somewhat more complicated by the existence of 3 mixing angles, 
two independent mass differences and many decoherence parameters. We follow an analogous 
method to that described above but now, instead of using the Pauli matrices as a basis, we 
choose the generators of SU(3). In order to derive analytical results, we choose the form of the 
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additional decoherence matrix so we may solve the differential equations straightforwardly but 
keep as many decoherence parameters as possible. Including standard oscillations, we therefore 
find the analogous sum of h + h' in Eq. (|38j) to be given by 
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(43) 



where Uji = ArrijjAE. Solving these equations and using Eq. 
bility from flavor p to q as 



gives the oscillation proba- 
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(44) 



where the C/'s denote the entries in the standard mixing matrix introduced in Eq. (j3J). If we 
again consider neutrinos travelling large distances from astrophysical objects, we can average 
the sin and cos terms to zero. If, for simplicity, we assume $ = 5, the probability reduces to 
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We can now derive equations describing the flavor composition at the detector in the same way 
as we did in the LV case: 
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[0.189e- 0.034]. 



(45) 
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C. CPTV and the LSND anomaly 

Having discussed how LV and CPTV may alter the phenomenology of neutrino oscilla- 
tions, we are now in a position to discuss these effects within the context of the anomalous 
LSND result. As we discussed in sectional the LSND experiment found a mass difference for 
anti-neutrinos which is incompatible with that of standard oscillations with three neutrinos. 
Assuming that these results are indeed correct, there are two ways to explain this result. The 
first is to assume that there are more than three neutrinos, with the additional neutrino(s) 
being sterile and not directly detectable. The second is to violate CPT invariance. 



1. Direct CPTV and the LSND anomaly 

In order to reconcile the LSND result with the rest of the neutrino oscillation data, it seems 
necessary to modify the neutrino sector to include different independent mass splittings for 
neutrinos and anti-neutrinos. In this way, we preserve the number of neutrinos but allow the 
mass differences to differ between the neutrino and anti-neutrino sectors. The present situation, 
taking in to account solar, atmospheric and KamLAND data, disfavours this scenario |57| . 
however, in both two and three neutrino models (Hil . IHgf , leaving room only for CPTV in four 
neutrino models where the mixing parameters may be different in neutrino and anti-neutrino 
sectors roC 



2. CPTV from quantum decoherence and the LSND anomaly 

A second possibility for explaining the LSND result without enlarging the neutrino sector 
is to consider quantum decoherence in the anti-neutrino sector only [61]. The oscillation prob- 
ability for three anti- neutrinos now takes the form of Eq. (J44j) . whilst neutrinos experience no 
quantum decoherence. Since decoherence parameters are present in the arguments of the sine 
and cosine terms, they alter the effective mass differences leading to an apparent difference 
between the measured mass differences in the neutrino and anti-neutrino sectors. It is therefore 
possible to reconcile the LSND results with other existing oscillation data 61]. However, this 
particular model fails to fit the spectral distortions observed in the KamLAND experiment 
|62j |. Having said that, the authors of Ref. [6]J chose only one set of quantum decoherence 
parameters and so there is still much scope for further investigation. 



V. SOURCES AND DETECTION OF HIGH-ENERGY NEUTRINOS 
A. Sources of high-energy neutrinos 

High-energy neutrinos are thought to be generated in a wide range of astrophysical sources. 
Such sources may include Active Galactic Nuclei (AGN) |63j|, Gamma- Ray Bursts (GRB) [64], 
microquasars |65[, supernova remnants, star clusters and X-ray binaries. Also, ultra-high pro- 
tons or nuclei travelling over cosmological distances can interact with the Cosmic Microwave 
Background (CMB) and/or Cosmic Infra- Red Background (CIRB), generating what is often 
called the cosmogenic neutrino flux (66[. 
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In any of these sources, there are basically two mechanisms by which high-energy neutrinos 
are generated. Firstly, Fermi accelerated protons (or charged nuclei) can collide with hadronic 
or photonic targets generating charged and neutral pions. These charged pions then decay, 
7r + — > /i + z/ M — > e + z/ e z/ e z/ M , 7T _ — > /i~z^t — > t~v e v e v^ generating electron and muon neutrinos and 
anti-neutrinos. Secondly, atomic nuclei which undergo Fermi acceleration can be disintegrated 
by interacting with infra-red photons surrounding their source. Neutrons broken off of such a 
nucleus can then decay, n — ■> p + e~v e , generating a flux of electron anti-neutrinos. 

The important thing to keep in mind regarding these two mechanisms for high-energy neu- 
trino generation is the quantity of neutrinos produced of various flavors. Neutrinos produced 
in charged pion decay follow the ratio: u e : : v T = 1/3 : 2/3 : 0, while those produced in 
neutron decay follow: v e : : v T = 1 : : 0. It is also important to note that pion decay can 
generate both neutrinos and anti-neutrinos, while neutron decay generates only anti-neutrinos. 

For the purposes of the measurement of flavor ratios, identified point sources of high-energy 
neutrinos are considerably more useful than diffuse fluxes. With such a source (or sum of 
sources), the distance the neutrinos have propagated will likely be known. Furthermore, the 
background from atmospheric neutrinos can be controlled by only considering events from one 
direction in the sky. In some sources which emit neutrinos only for short lengths of time, GRB 
and AGN flares for example, the background can be further reduced by only considering events 
in particular time windows. For these reasons, GRB and AGN are likely to be among the most 
useful source of high-energy neutrinos for the purposes of flavor identification although bright 
and nearby (galactic) sources, if present, could also be very useful. 

In addition to these theoretical arguments, there is some limited experimental evidence that 
might suggest the existance of bright point sources of high-energy neutrinos. Firstly, it has been 
argued recently that anisotropics observed in the cosmic ray spectrum at EeV energies is the 
result of neutrons propagating from galactic sources. If this is the case, then large fluxes of high- 
energy (anti-)neutrinos will also be generated |67J. Secondly, the AMANDA-II experiment has 
recently reported the detection of two neutrinos coincident with TeV flares seen by the Whipple 
gamma-ray telescope from the blazar (AGN) 1ES 1959+650 [68]. These events were not found 
in a blind analysis, however, so their statistical significance cannot be determined. If these 
neutrinos are the product of this TeV blazar, it would suggest the existance of very bright 
point sources of high-energy neutrinos. 

More exotic processes which do not fall into these two descriptions may also be capable of 
generating high-energy cosmic neutrinos. Such possibilities include annihilating or decaying 
dark matter or topological defects Hawking radiating primordial black holes [fjj, or the 
interactions of ultra-high energy neutrinos with the cosmic neutrino background via the Z- 
burst mechanism [7]]]. For a review of sources of high-energy neutrinos and other aspects of 
high-energy neutrino astronomy, see Ref. [72l ]. 

B. High-energy neutrino detection 

Once such neutrinos are generated and propagate to Earth, they can be detected in one of 
several ways. Neutral current interactions of neutrinos of all flavors with nucleons generates 
hadronic showers which can be observed. Charged current interactions of electron and muon 
neutrinos generate, in addition to hadronic showers, potentially observable electromagnetic 
showers and muons, respectively. The tau leptons generated in the charged current interactions 
of tau neutrinos can produce a class of events unique to tau neutrinos: double bangs and 
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lollipops. 

Many of the experimental techniques being developed and deployed are only capable of 
detecting showers generated in high-energy neutrino interactions. Observing shower events 
alone will not enable the flavor ratios of a flux of cosmic neutrinos to be identified, however. 
For this reason we focus on experiments capable of observing high-energy neutrinos in the 
form of showers, muon tracks and tau-unique events. In particular, we focus on next gen- 
eration, kilometer-scale, optical Cerenkov detectors. These include IceCube, currently under 
construction at the South Pole, and possibly a future kilometer-scale neutrino telescope in the 
Mediterranean Sea, sometimes called KM3. 



1. Shower events 

All high-energy neutrino interactions produce an electromagnetic and/or hadronic shower. 
The probability of detecting a hadronic shower produced in a neutral current interaction as a 
neutrino travels through the effective area of the detector is given by 

-P^showcr = pN A L / -i-dy, (46) 
Je$/e v d y 

where p is the target nucleon density, N A is Avogadro's number, L fa 1 km is the leng th of the 
detector, da/dy is the differential neutrino-nucleon neutral current cross section [73j|, y is the 
fraction of energy which is transferred from the neutrino (and therefore the fraction of energy 
which goes into the shower) and E^ fa 3 TeV is the threshold energy for the experiment 
detecting a shower event. 

In the charged current interactions of electron neutrinos, all of the neutrino's energy goes 
into a combination of electromagnetic and hadronic showers. In this case, the probability of 
detecting a shower reduces to 

^shower = pN A La, (47) 

if E v > E^, and zero otherwise, a in this expression is the total charged current neutrino- 
nucleon cross section 

0- 



In these expressions, we have made the simplifying assumption that only showers generated 
inside of the detector volume can be detected. This is not always true, particularly at very high 
energies. Very energetic showers which are initiated outside of the experiment's instrumented 
volume can expand into the experiment. Another way of saying this is that the effective volume 
of such an experiment is generally larger than its instrumented volume at very high energies. 



2. Muon events 

Charged current interactions of high-energy muon neutrinos produce muons which can travel 
through the medium of the experiment (ice or water) and potentially into the detector volume. 
The energy loss rate of such a muon is given by 

dE 
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where the parameters in ice or water are given by a « 2.0 MeV cm 2 /g and (3 « 4.2 x 10" 



cm 2 /g |74J]. The distance a muon travels before its energy drops below the threshold, E^, is 



given then by 



R„ 



a + (3Ef 



(49) 



fj, j 

This quantity if often referred to as the muon range. E^ v ~ 50-100 GeV are typical for high- 
energy neutrino telescopes. The muon range can extend for many kilometers for very high 
energy muons, dramatically increasing the number of muon tracks that are observed. 

The probability of detecting a muon produced in a charged current interaction is given by 



P N A aR 



in 



(50) 



where, here, a is the total charged current neutrino-nucleon cross section 73]. In cases with 
the muon neutrino comes from above the ice or water, the value of used in this equation 
should not exceed the depth of the experiment. Similarly, if the neutrino comes from below the 
detector, the rock or other material below the ice or water should be accounted for. 



3. Events involving tau neutrinos 

For tau neutrinos with energies less than ~PeV, their interactions generate only shower 
events, as any tau leptons generated decay before they can be identified. Higher energy tau 
neutrinos, on the other hand, can be identified by the combined signatures of tau tracks and 
showers. 

Double bang 7^, 76[ events are produced with a tau neutrino interacts via charged current 



inside of the detector, producing a tau lepton which travels across the detector to decay, still 
inside of the detector volume. The hadronic shower produced in the initial interaction consti- 
tutes the first "bang" while the second shower is generated in the tau's decay. A tau lepton 
can travel a distance 

R T = — — - ~ r — 50 meters (51) 
m T 1 PeV 

before decaying. To resolve two showers in a high-energy neutrino telescope, they must be 
separated by a distance of roughly 100-400 meters. Furthermore, they cannot be separated by 
more than about 1000 meters and both be within the detector volume. Therefore double bang 
events are most useful in the energy range between a few PeV and one hundred PeV. 

At higher energies, lollipop events become very useful [76j]. A lollipop event is observed when 
the first shower of a double bang event occurs outside of the detector (without being observed) 
with the tau track extending into the detector and decaying. You might also think that an 
observation of the first shower with a tau track would be a useful signature, but muons and 
taus produced in ordinary hadronic or electromagnetic showers could mimic such an event. 

The probability of observing a double bang event from a given incident tau neutrino is given 

by 

Pdb « P N A f dy^- [ L dx { -^±e-^, (52) 
Jo dy A min Rt 



where as the probability for a lollipop event is 

f dy 

o dy 



> LP « P N A {L - x min ) I dy^e~ x ^. (53) 
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In each of these expressions, x m - m ~ 100 — 400 meters, is the minimum shower seperation 
required to seperate the showers. 

For a detailed description of both double bang and lollipop tau events, see Ref . [76| . 



4- Neutrino events in cosmic ray experiments 

Although in this article we focus on next generation optical Cerenkov neutrino telescopes, it 
is interesting to point out that ultra-high energy cosmic ray experiments may also have a limited 
ability to resolve neutrino flavors. These experiments can detect ultra-high energy neutrinos in 
essentially two ways. First, quasi-horizontal neutrinos can penetrate deeply into the atmosphere 
before interacting, producing showers which are distinguishable from those initiated by cosmic 
rays. Second, tau neutrinos which skim the Earth can, through charged current interactions, 
produce tau leptons which escape the Earth before decaying and generating a shower |77| . 
While the former signature can be produced by neutrinos of all three flavors, the latter are 
uniquely generated by tau neutrinos. This, in principle, could be used to measure the fraction 
of the ultra-high energy neutrino flux which consists of tau neutrinos. 



VI. SIGNATURES OF CPT AND LORENTZ VIOLATION 

There are two primary methods of probing physics beyond the Standard Model with high- 
energy neutrinos. First, neutrino-nucleon interactions can be observed measuring the cross 
section [78| and other characteristics in hope of identifying new interactions: microscopic black 
hole or P-brane production [7^|. processes resulting from low-scale gravity string effects 
jU or electoweak instantons |82|. Second, the ratios of neutrino flavors can be measured, 
potentially identifying the effects of neutrino decay [jjj . pseudo-Dirac states [§i| or quantum 
decoherence |52j ]. In this article, we focus on using this latter technique to constrain or discover 
the effects of CPT and Lorentz violation. 



A. Flavor ratio predictions 

To assess the prospects of identifying the violation of CPT and Lorentz invariance in the high- 
energy neutrino sector, we must first determine the flavor ratios predicted in various scenarios. 
First of all, we consider the case with only known physics, and no CPT or Lorentz violation. For 
an initial set of ratios corresponding to neutrinos from pion decay, v e : : v T = 1/3 : 2/3 : 0, 
after oscillations over a long distance (which will always be the case for high-energy neutrino 
astronomy), these ratios become u e : : v T « 0.36 : 0.33 : 0.30. (Anti)-neutrinos coming from 
neutron decay, on the other hand, are generated in the ratio u e : : v T = 1 : : 0, which 
oscillates to v e : : u T m 0.56 : 0.26 : 0.18. 

These ratios can be changed dramatically if the effects of CPT or Lorentz violation are 
present. In the case of LV, above the energy threshold for such effects, the neutrino flavor 
ratios from pion decay are modified as v e : : v T = 1/3 : 2/3 : — > 0.42 : 0.57 : 0.013. The 
flavor ratios from neutron decay are modified as v e : : v T — 1 : : — ► 0.70 : 0.27 : 0.027. In 
the case of quantum decoherence, on the other hand, all ratios shift toward v e : : v T — 1/3 : 
1/3 : 1/3, regardless of their source (see figure |2J). 
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B. Lorentz Violation 



There are two potentially identifiable features in the neutrino flavor ratios predicted as 
a result of LV. First, in the case of neutrinos generated via pion decay, a particularly large 
fraction of these neutrinos will be of muon flavor. Secondly, regardless of whether the neutrinos 
are generated in pion or neutron decay, only a very small fraction will appear with tau flavor. 

1. A large muon neutrino fraction 

If the effects of LV are considered in neutrinos generated via pion decay, the result can be a 
neutrino flux which is of nearly 60% muon flavor. A large muon to shower ratio could, therefore, 
be seen as a signature of LV. The ability of high-energy neutrino telescopes to measure the ratio 
of muon to shower events and corresponding flavor ratio is discussed in detail in Ref. |76^ . 

2. A tau neutrino deficit 

Looking for a deficit of high-energy cosmic tau neutrinos may also be useful in identifying the 
effects of LV. Events that are identifiable as tau neutrinos (double bang and lollipop events) are 
somewhat rare, however. Assuming a spectrum proportional to E~ 2 , a flux of E 2 T dN Ur j dE VT = 
10 -7 GeV cm -2 s~ x would yield only about 0.35 tau-unique events per year in a cubic kilometer 
experiment. Over several years of observation, observing a deficit of tau neutrinos may indeed 
become possible if bright point sources of PeV-EeV neutrinos are discovered. 

The absence of tau neutrino induced Earth-skimming showers at cosmic ray experiments, 
such as Auger, could also be an anticipated signature of LV effects. As experiments such as 
Auger most efficiently detect showers with energies above ~ 10 8 GeV, or so, they are particularly 
well suited for testing the effects of LV. 

C. Quantum Decoherence 

The signature of quantum decoherence in cosmic neutrinos, is the presence of an equal 
fraction of neutrinos of each flavor. Since sources which produce neutrinos through pion decay 
generate nearly this ratio (after oscillations) without the effect of quantum decoherence, they 
are not very useful in probing for these effects. Sources of high-energy (anti-) neutrinos produced 
through neutron decay, on the other hand, can be used to potentially identify these effects |H3 |. 

In figure EJ we show the effects of quantum decoherence on the events observed in a next 
generation high-energy neutrino telescope. In the left frame, we consider a model with a 
5 oc E 2 dependence, normalized such that quantum decoherence sets in at the 10 TeV scale 
(5 = (E/10, 000 GeV) 2 /L). In the right frame, the source is sufficiently distant that the effects 
of quantum decohence have set in fully [y e : : v T = 1/3 : 1/3 : 1/3). In either case, the 
number of muon tracks detected are increased and the number of showers detected decreases. 
Since the normalization of this flux will likely be unknown, it is the ratio of these event types 
that is of the most interested to us. 

In figure H3 we plot the total rate of muon events and the rate of contained muon events (i.e. 
with a vertex contained within the detector volume) seperately. We do this to illustrate that 
energy dependent effects, such as those shown in the left frame of this figure, are more clearly 
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FIG. 2: The effect of quantum decoherence on the ratios of neutrino flavors. In the left frame, the 
initial ratios are set to those values found from pion decay {y e : Va : v T = 1/3 : 2/3 : 0) while in 
the rigth frame the ratios from neutron decay (y e : Va : v T = 1 : : 0) are shown. In both cases, 
we have considered a model with 5L = (£7/10, 000 GeV) 2 . In the case of neutrinos generated in pion 
decay, quantum decoherence has little impact on the flavor ratios (left frame). The flavor ratios of 
(anti-)neutrinos generated in neutron decay, on the other hand, can be dramatically effected (right 
frame) . 



identified by using only contained events. The statistics are considerably better when using all 
of the events, however. Depending on the flux of high-energy neutrinos present, a strategy to 
use these different types of events most effectively will have to be devised. 

The prospects for detecting the effects of quantum decoherence or other signatures of CPT 
or Lorentz violation ultimately depend on the variety of high-energy neutrino sources which 
exist in nature. A detailed study of IceCube's sensitivity to decoherence effects from observing 
anti- neutrinos from the Cygnus spiral arm is currently underway j85[. This is likely to be one 
of the most useful sources for constraining the effects of decoherence. 



VII. CONCLUSIONS 

High-energy neutrino astronomy provides an opportunity to observe particles at extremely 
high energies and over extremely long baselines. Both of these characteristics make such ex- 
periments particularly adept at testing for the effects of CPT and Lorentz violation. 

In this article, we have discussed in detail the effects that CPT and Lorentz violation can 
have on the flavors of high-energy cosmic neutrinos observed at Earth. After discussing the 
theoretical basis and motivations for such effects, we calculated the ratios of neutrino flavors 
predicted to be observed from astrophysical sources of high-energy neutrinos in various CPT 
and Lorentz violating scenarios. 

The effects of Lorentz violation may potentially be detected or constrained by the observa- 
tion, or lack thereof, of anomalously large fractions of the cosmic neutrino spectrum consisting 
of muon neutrinos, accompanied by very few tau neutrinos. This deviation from the standard 
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FIG. 3: The effect of quantum decoherence on the events observed in a next generation high-energy 
neutrino telescope, such as IceCube or KM3. Solid lines are the predicted distribution of events with no 
effects of quantum decoherence, while the dashed lines represent the cases of 5L = (22/10, 000 GeV) 2 
(left frame) and full decoherence (right frame). The (from top to bottom on the left side of the 
frames) blue, red and black lines show the distributions of muon events, electromagnetic/hadronic 
shower events, and all muon events with a vertex contained in the detector volume, respectively. The 
flavor ratios used are for a source producing (anti-)neutrinos via neutron decay. The rates have been 
normalized to a total neutrino flux of E„dN v / dE u = 10~ 7 GeV cm -2 s _1 . 



neutrino oscillation phenomenology will occur above a model-dependent energy threshold, and 
thus experiments capable of detecting extremely high energy neutrinos are particularly useful 
in such measurements. 

The effects of quantum decoherence may also be observable in cosmic neutrinos. In particu- 
lar, anti-neutrinos generated in the decays of neutrons produced in the photo-disintegration of 
ultra-high energy cosmic nuclei will have their flavor ratios significantly modified if quantum 
decoherence effects are significant. The very long baselines over which such neutrinos travel be- 
fore reaching Earth provide an opportunity to test for these effects with much greater precision 
than other techniques can achieve. 
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